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Abstract 

In the framework of quantum electrodynamics (QED) in external 
potentials, we introduce a method to compute the time-dependence of 
the expectation value of the current density for time-dependent homo¬ 
geneous external electric fields. We apply it to the so-called Sauter 
pulse. For late times, our results agree with the asymptotic value due 
to electron-positron pair production. We correct, and compare to, a 
general expression derived by Berber for the linearization in the ex¬ 
ternal field. Based on the properties of the current density, we argue 
that the appearance of enhanced quasi-particle densities at intermedi¬ 
ate times in slowly varying sub-critical potentials is generic. Also an 
alternative approach, which circumvents these difficulties, is sketched. 

Classification (PACS): 12.20.Ds; 11.40.-q 
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1 Introduction 

A crucial prediction of QED in external potentials is the production of 
electron-positron pairs in sufficiently strong electric fields III Eli, the so- 
called Schwinger effect. This is a non-perturbative effect which is exponen¬ 
tially suppressed for field strengths below the critical field strength Ec = 
w?/\e\ ~ 10^®V/cm. Hence, it seems unrealistic to test it with static fields. 
More promising are high intensity lasers, cf. the reviews and references 
therein. However, even the projected high intensity laser facilities may not 
reach this non-perturbative regime. One possibility to overcome this is to 
use complicated field configurations, such as in the dynamically assisted 
Schwinger effect [6]. Often, such configurations can only be treated numer¬ 
ically. One framework for this is the quantum kinetic equation (QKE) [7j, 
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where one defines quasi-particle number densities at intermediate times and 
obtains the asymptotic particle number by solving an ordinary differential 
equation for their time evolution^ It was noticed [H], that the quasi-particle 
numbers, as defined in this framework, can be much larger at intermedi¬ 
ate times than at asymptotically late times (for parameters typical for an 
optical laser setting, the proportionality was 10^^). This is problematic in 
numerical studies, as the errors that are made at intermediate times, where 
the number densities are high, are propagated to asymptotic times, where 
they can dominate the actual number densities. 

Now the concept of a particle is ambiguous at intermediate times, where 
the external potential is non-vanishing and varying; consequently, one speaks 
of quasi-particles^ One may thus wonder how generic the occurrence of en¬ 
hanced quasi-particle densities at intermediate times is. Or more generally: 
consider any formalism in which particle creation is studied by (numeri¬ 
cally) solving an evolution equation for some dynamical quantity]^ Can one 
make any generic statements about the intermediate time behavior? One 
possibility to investigate this is to consider a quantity that is unique and 
well-defined also at intermediate times. The prime candidate for this is the 
current densitj]^ 

j^^ = ( 1 ) 

It is an observable, and thus, at least in principle, measurable. Such a 
measurement would proceed via back-reaction, i.e., by the changes in the 
electromagnetic field which are sourced by this current. These back-reaction 
effects can be used to constrain the generation of primordial electromagnetic 
fields during inflation m- However, in the present work we are more inter¬ 
ested in generic features of this observable and whether these may help to 
better understand the intermediate time behavior of the dynamical variable 
in the time evolution equation in which the system is studied. In particu¬ 
lar, we will argue that the occurrence of enhanced quasi-particle numbers 
at intermediate times is generic, and will also indicate a new scheme which 
avoids these difficulties. 

The current density 0 is unique in the following sense: For a locally 
gauge covariant definition of the current density, the only ambiguity consists 

'^To be precise, in order to avoid an integro-differential equation, one either has to 
supplement the quasi-particle number density with other variables, as in 0 below, or use 
the underlying Bogoliubov coefficients, cf. [8] for a detailed discussion. 

^For a recent discussion of these ambiguities in the context of scalar QED in external 
potentials, we refer to [10]. In the context of quantum held theory on curved space-times, 
these issues were discussed for example in m m- 

^Another such framework is the Wigner function formalism introduced in |13| . cf. mi, 
for example. Compared to the quantum kinetic equation it has the advantage of being 
in principle applicable also when the external held is more general than a homogeneous 
electric held. 

"^The analog of this for quantum held theory in curved space-times would be the stress- 
energy tensor. 
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in adding multiples of the external current density uni, corresponding to a 
charge renormalization. It can be fixed by requiring that for static external 
potentials the total vacuum polarization charge vanishes at linear order in 
the external potential [IZIE], i.e., the external charges are not renormalized 
at this order of approximation. 

The expectation value of the current density Q, henceforth also called 
the vacuum current density, is a crucial ingredient in the study of back- 
reaction effects. In the QKE framework, an expression for the vacuum cur¬ 
rent density was given in m, cf. also m- In fact, the analysis of this 
expression will be crucial for the explanation of the enhanced quasi-particle 
numbers in slowly varying sub-critical potentials reported in [9]. However, 
the renormalization prescription used there is not explicitly local, so it is 
a priori not clear whether it conforms to the requirement of local gauge 
covariancej^ It it thus desirable to check these results with independent 
methods. 

Based on Dirac’s method of point-splitting w.r.t. the singular part of 
the two-point function]^ we give, in Section]^ a general expression for the 
vacuum current in a homogeneous time-dependent external electric field. We 
show that it is equivalent to the expression derived in the QKE framework. 
We also evaluate it for the concrete case of the so-called Sauter pulse, 

E{t) =- K —, A{t) =-^tanhrct, (2) 

cosh wt w 

both for super- and sub-critical peak field strengths Eq. This external field 
was studied extensively in the literature, with an emphasis on the asymptotic 
behavior, i.e., the pair production probability, which was first computed by 
Narozhnyi and Nikishov [23]. The intermediate time behavior of the quasi¬ 
particle number in this system was studied in [U |25| using the quantum 
kinetic equation |7]. Our results for the Sauter pulse agree well with the 
asymptotic pair production probabilities derived in |24) . 

We also recall, in Section]^ that there is a general expression for the vac¬ 
uum current density at linear order in the external potential, due to Serber 
|26j . However, his derivation contains some mistakes, leading to an incor¬ 
rect result, which we correct here. For sub-critical peak field strengths, the 
application of this general expression to the Sauter pulse Q leads to a nice 
agreement with the results obtained by our new method. One interesting 

®We refer to [20] for an elementary example showing that in the presence of external 
potentials a seemingly sensible but not explicitly local renormalization prescription can 
give wrong results. 

®This is the analog of the covariant point-splitting renormalization of the stress-energy 
tensor in the context of quantum field theory in curved space-times m, m\, cf. also [23] 
for a recent review. 

^In the former reference, no physical interpretation of the quasi-particle number at 
intermediate times was attempted, noting that an interpretation as real particles is only 
possible at asymptotic times. 
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feature of Berber’s general expression is that, at the linear order, the vacuum 
current only depends on the external cnrrent, not on the field strength. 

We remark that the study of generic properties of the current density 
does not seem to have attracted a lot of attention in the literature. Usu¬ 
ally, the current density is considered in the context of back-reaction, with 
an emphasis on directly including these effects in the quantum evolution 
equation For a work where the vacuum current was computed 

independently of back-reaction, we refer to |28j . 

Finally, in Section we turn to the discussion of the enhanced quasi¬ 
particle numbers reported in [9] for slowly varying sub-critical potentials. 
For the adiabatic basis nsually considered, this turns out to be a straight¬ 
forward consequence of the expression for the vacuum current density in 
terms of the qnasi-particle density. We argue that this is indeed a generic 
feature. Also a new framework, based on Dirac’s point splitting method, is 
sketched, in which these problems should be absent, and which is in principle 
also applicable for generic external fields. 

2 The current density 

For the computation of the vacuum current density we use Dirac’s method 
|29j of point-splitting w.r.t. the singular part of the two-point function 
('0(x)'0(y)), the so-called Hadamard parametrix H{x,y), cf. [20] for a de¬ 
tailed discussion. This method is conceptually appealing as it is explicitly 
local and gauge covariant. The practical difficulty is that the Hadamard 
parametrix is known in position space, while the two-point function is typi¬ 
cally given by an integral over modes. The subtraction of the singular part 
H from this mode integral is in general quite challenging. In |20| a solntion 
to this problem for the computation of the vacuum polarization in static 
external potentials was proposed. The idea is to consider the limit of coin¬ 
ciding points from the time direction and to rewrite + se*^) (up 

to a potential logarithmic singularity in s) in the form of an integral over 
(shifted) vacunm mode densities. The snbtraction of these singular parts 
can thus be performed inside the mode integral, which renders it much less 
divergent and numerically more tractable. 

This method can be adapted to the case of a homogeneous external 
electric field along the z-axis, in the gange A^^{x) = 5^A{x^). It is then 
advantageous to perform the limit of coinciding points from the z-direction. 
An analogous rewriting of the Hadamard parametrix yields the following 


4 


expression for the vacuum current; 


= e lim 

z —>-0 


dp 3 {fiP3,x°) - fo{p3 - eA{x^), x°)) 

+^ej3(x°) (log.VA^-1) 


( 3 ) 


Here / is defined as follows. Consider the properly normalized mode solu¬ 
tions of the Dirac equation 

— m)ilj = 0 


in the presence of the external potential, which coincide with the posi¬ 
tive/negative energy modes at asymptotically early times (we refer to Ap¬ 
pendix for a precise definition of the required normalization). Here s is a 
spin label, p a three-momentum, and x = {x^,x). Then f{p 3 ,x^) is defined 
as 

f{p3,x°) = '^ [ dV fpAx^)7^fp,s{x°), 

where p±_ = {pi,P 2 ) is the perpendicular momentum. The function /o in 
Q is defined analogously, but for the standard mode solutions fp]t in the 
absence of external potentials. The corresponding term in (§ can thus be 
seen as the subtraction of the vacuum contribution, with an appropriate shift 
in the momentumj^ The third term in Q cancels a remaining logarithmic 
divergence in the presence of an external current A = d\F^^. As discussed 
in the introduction, we fix the length scale A such that external charges are 
not renormalized at the linear level, i.e.. 


A = 2/(e'’'m), 


( 4 ) 


cf. [2], with 7 the Euler-Mascheroni constant. In fact, using the relation 
([T^, Q may be rewritten as 


(fix)) 

where 




o,+ 

7r,s 


(27r)3 2Aj ’ 


^1/2—Pi/25 (5) 

vra =P3 - eA(x°), 

2 2 I 2 I 2 

Sj_ = m -hpi +P2, 

ujp = nA + \pA 


®One may wonder whether there is a gauge dependence in due to the occurrence 
of A{x^). However, a shift A ^ A + a, has to be accompanied by the replacement 
fiP 3 ,x°) f{p 3 -ea,x°) in as the new solution to the equation of motion ( |12| for 
the mode functions is given by fA = /j(-eae3 s- With the substitution qa = ps — ea, this 
results in a phase shift of the integral in which is irrelevant in the limit 2 —>■ 0. 
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There is also an expression for the vacuum current density in the frame¬ 
work of the QKE, 


(fit)) = -2e 


d^p 

(27r)^ 


— n{p) + 

UJt^ 


CJ71 




eE{t)e‘j 


( 6 ) 


due to [l8], cf. also |8]. Here n{p,t) is the number density of quasi-particles 
of momentum p at time t, fulhlling mu 


h{p,t) = Q{p,t)k{p,t), 

k{p,t) = Q{p,t){l - n{p,t)) - 2uj{p,t)£{p,t), (7) 

e{p,t) = 2uj{p,t)k{p,t), 


with 


Q{p, t) 


eE{t)s± 


and the initial conditions n{p,to) = k{p,to) = £{p,to) = 0 for some time to 
before the potential is switched on. 

The hrst term in the integrand in Q gives rise to the so-called conduction 
current, due to the movement of the quasi-particles already created. The 
last two terms yield the polarization current, due to the creation of quasi¬ 
particle pairs. The last term is a counterterm, canceling a divergence in the 
polarization current. This will be crucial in our discussion of the enhanced 
quasi-particle numbers at intermediate times. Note that the expression Q 
for the vacuum current is the correct expression for quasi-particle densities as 
defined in the QKE framework. For other definitions of quasi-particles, the 
expression has to be modihed (but note the discussion in Section below). 

As discussed in Appendix the equations (© and are equivalent, 
the precise relation being given by @. 

The mode solutions for the Sauter pulse (§ are explicitly known 
|24| . so the expression Q can be evaluated numerically rather straightfor¬ 
wardly]^ For super-critical peak field strengths Eq > Ec, we find that the 
vacuum current density at intermediate times increases quasi-monotonically 
to its asymptotic value, even for very short pulses, cf. Fig. Furthermore, 
the slope of the vacuum current density is largest approximately where the 
field strength has its maximum, as one might expect. Note that also the 
asymptotic value is shown, as computed by integrating the pair production 
probability density P{p) of Narozhnyi and Nikishov |24) with the group ve¬ 
locity v{p). The vacuum current density nicely converges to this asymptotic 


®This is especially the case in the absence of a hierarchy of scales, l.e., when the 
“frequency” w is of the same order of magnitude as the mass m. For simplicity, we 
restrict to that case. Investigations with longer pulse durations show that the results 
are qualitatively similar, with the exception of the damping of the oscillating tail in the 
sub-critical case, cf. Fig. and the discussion at the end of Section 
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Figure 1: The vacuum current density (blue discs) as given by Q for the 
super-critical peak field strength Eq = lOEc and m = w = 1. The red 
line indicates the asymptotic value. Also shown are the conduction (violet 
squares) and polarization (yellow diamonds) current densities, cf. Q. For 
comparison, also the external cnrrent density, divided by 100, is shown (blue 
dashed curve). All quantities are expressed in natural Compton units. 


value. Note that, as indicated in the figure, the vacuum current density 
is two orders of magnitnde smaller than the external one. Hence, the ne- 
glection of back-reaction effects seems justified, at least on the time scales 
considered herej^ Also shown in the figure are the conduction and polar¬ 
ization currents as computed by Q using the quasi-particle densities n{q, t) 
for the Sauter pulse [8]. 

For sub-critical peak field strengths Eq < Ec, we hnd that the vacuum 
current density at intermediate times can be much larger than at asymptot¬ 
ically late times, cf. Fig. where the asymptotic current density is about 
three orders of magnitudes smaller than at intermediate times. Furthermore, 
for very short pulses, it oscillates after the passing of the pulse with a period 
of about half the Compton time of the electron. Now the vacuum current 
as defined in the QKE framework is nearly entirely due to polarization (so 
that the decomposition into the conduction and polarization current is not 
shown, the polarization current being indistinguishable from the full result). 

[TB] a seemingly different conclusion was obtained. However, the electric charge 
was there set to = 4 instead of the value = 47r/137 used here. 
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Figure 2: The vacuum current density (blue discs) as given by Q for the 
sub-critical peak field strength Eq = Q.lEc and m = w = 1. Now the red line 
is the prediction from Berber’s equation Q. The external current density, 
divided by 1000, is given by the dashed bine line. 

3 Linearization in the external field 

One of the first applications of QED in external potentials was Uehling’s 
calcnlation of vacnum polarization in static potentials, at linear order in 
the external potential |30j . A generalization for the time-dependent case 
is due to Berber [26]. However, his derivation is not entirely correct, cf. 
Appendix]^ A corrected version is 

{j^{x)) = j d^x' K{x — (8) 

where J is the external current density and K vanishes outside of the forward 
light cone and is given by 

^ poo 

K{x) = - -j= / d'i/’ cos^V’ / dk Ji {2kV^/cos (9) 

Svr^vx^ Jo Jm 

inside, with Ji being a Bessel function. Interestingly, ([^ only depends on the 
external current, not on the field strength. Heuristically, in time-dependent 
situations, the vacnum current Q can be thought of as the current dne to 
a rearrangement of the Uehling vacuum polarization clouds accompanying 
the charges that are responsible for the external current. In this sense, it 
is indeed a polarization cnrrent, as in the terminology of [T8|, discussed 
above. However, it may be misleading to attribute it to pair creation, as 
usually done. An expression similar to ([^ was found by Bchwinger m, 
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Figure 3: Berber’s kernel K{x), cf. multiplied by Distances are 
measured natural Compton units. 


who, however, assumed an external potential vanishing in the past and the 
future, a condition not fulfilled in the present case. 

As shown in Figure]^ the expression Q describes the time-dependence 
of the vacuum current in the sub-critical case very well. We note that 
in the example shown in Figure the adiabaticity parameter mu}/\eE\ is 
larger than one, indicating that the (asymptotic) pair production probabil¬ 
ity is dominated by multi-photon processes [32]. Nevertheless, the vacuum 
current density at intermediate times is very well described by the linear 
approximation. We also note that the asymptotic value of Berber’s expres¬ 
sion Q for the vacuum current density vanishes. This follows from the 
Riemann-Lebesgue lemma and the integrability of the Fourier transform of 
K at its singularity = 4m^. 

The oscillatory behavior of the current shown in Fig.[^is a manifestation 
of the oscillatory nature of K, cf. Fig. Clearly, for external currents that 
vary on time-scales much larger than the Compton time of the electron, 
these oscillations are strongly damped. Following Bchwinger m , one may 
then write the vacuum current density as 

oo 

n=l 

for certain constants Cn (which can also easily be obtained by expanding the 
function x(/c) of Berber |26| around k = 0). 
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The enhanced quasi-particle densities at inter¬ 
mediate times 


Let us now turn to the discussion of the enhanced quasi-particle densities 
described in [9]. There, a slowly varying sub-critical electric field E = 
Eosinwt is considered, which is turned on and off at some finite times. At 
first sight, one might think that the explanation lies in the fact that for sub- 
critical electric fields the vacuum current density at intermediate times is 
dominated by the linearization Q, which is known to vanish at asymptotic 
times. A first indication that this reasoning is flawed is the fact that one 
may well have non-vanishing particle densities at asymptotic times in spite 
of a vanishing vacuum current. As we will see, a satisfactory explanation of 
the enhanced quasi-particle densities at intermediate times has to take the 
logarithmic divergence of the current into account. 

We recall that the third term in Q is a counterterm which is necessary to 
cancel a logarithmic divergence in the polarization current, i.e., the integral 
over the second term in Q . It follows that for large momentum p the rate 
of the change of the quasi-particle density is given bj 
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h{p,t) ~ e^E{t)E{t) 


4w(p, tf 


For slowly varying sub-critical external fields, one expects that this yields 
the main contribution to the rate of change of the total quasi-particle den¬ 
sity. Performing the momentum integral over the fraction on the r.h.s. and 
integrating in time, this implies 




( 10 ) 


Up to a factor ^ ~ 0.93, which seems to have been interpreted as 1, 
this agrees with the heuristic description of the numerical results found in 
|9j for a sinusoidal time-dependence of the background electric field. The 
expectation (10) can also be verified for slow sub-critical Sauter pulses ([^, 
using the quasi-particle number density from |8], cf. Figure]^ In this sense, 
the enhanced quasi-particle densities reported in [9] are primarily due to 
the fact that the dynamical variables chosen in the QKE framework are 
such that a supplementary counterterm is necessary in the expression for 
the vacuum current density, i.e., that h is much larger than it would need 
to be in order to yield the correct vacuum current density. 

A strategy to circumvent the problem of greatly enhanced particle num¬ 
ber densities would thus be to find dynamical variables for which no further 
renormalization is necessary in the computation of the current density. One 
possibility would be to consider better adapted modes for the definition of 


Actually, in the counterterm in (§ is derived from this asymptotic form of h. 
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Figure 4: The quasi-particle number density as defined in the framework 
of the quantum kinetic equation (blue discs), compared to the expectation 
(10), for the Sauter pulse ([^ with Eq = O.lEc, w = 0.1, m = 1. 


quasi-particle densities, perhaps similar to those discussed in m for the 
charged scalar field. However, it seems unlikely that this is possible: we 
know that the /3 function of QED does not vanish, so there must be a 
logarithmic divergence in the current, proportional to the external current 
J = —E (a change j —>• J -b AJ corresponds to a charge renormalization). 
Hence, for any definition of quasi-particle densities, one expects a logarith¬ 
mic divergence proportional to E in the corresponding expression for the 
vacuum current. 

Instead, one might proceed as follows: The essence of Dirac’s procedure is 
to subtract from the two-point function (V'(x)'i/’(y)) its singular part y), 
the Hadamard parametrix. The result R{x, y) is smooth, so that the limit 
of coinciding points may be performed, for example in the computation of 
the expectation value of the current density or the stress-energy tensor. The 
idea would be to solve the differential equation for the smooth remainder 
R0i.e., 


DR{x,x') = DH{x,x'), 


D'^R{x, x') = D'^H{x, x'), (11) 


where D'^ is the transpose oi D = i'f — m (acting on a co-spinor), applied 
on the primed variable. Here we used that the two-point function is a 
bi-solution to the Dirac equation. In contrast, the parametrix is not a bi¬ 
solution, but there are smooth remainders which appear as source terms 
on the right hand side. These remainders are given [33] as a series in the 


may be advantageous to restrict to coinciding times and to include a compensating 
parallel transport in the definition of R, similarly to the Wigner function framework [13] . 
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Hadamard coefficients Vfc(x, x') (the Lorentzian analogues of the heat kernel 
coefficients) for the wave operator 

P = + f 


Instead of a linear differential equation for the Bogoliubov coefficients, one 
thus has to solve an inhomogeneous one for the remainder R. The renor¬ 
malization ambiguity of the current is now hidden in the initial condition 
for R (choosing P = 0 in the past amounts to Heisenberg’s renormaliza¬ 
tion condition Q)- The vacuum current density {j^{x)) is simply given by 
tr 7 ^P(x,x), i.e., no further renormalization is necessary. For the solution 
R one would thus expect that the magnitude at intermediate times is of the 
same order as asymptotically. Once the background field is switched off, 
and hence H equals the vacuum two-point function, one can read off the 
particle content directly from the Fourier transform of R. 

A nice feature of this scheme is that, due to the initial condition P = 0, 
all pair creation effects are due to the inhomogeneous part on the r.h.s. of 
0. Hence, it may be possible to make statements about pair creation 
without actually solving the differential equations (11) A further advan¬ 


tage of the framework is that it is in principle applicable to generic external 
potentials, in contrast to the QKE, which requires a homogeneous, purely 
electric field. 

On the other hand, there are obvious difficulties: the construction of H 
involves a series 

^ Vk{x,x'){x - x'f^^~^\ 

k 


which is only guaranteed to converge for analytic external fields. In practice, 
one would probably truncate the series, so that P is not smooth (but for 
arbitrarily large m). This potentially leads to problems at large distances, 
so it might be necessary to include an infrared cut-off in the definition of 
H. A further difficulty lies in the fact that the Hadamard coefficients 14 
are not known explicitly to all orders (however, for very symmetric external 
fields it could be possible to fully determine them). In view of the potential 
benefits, we consider it worthwhile to try to overcome these difficulties. 

We remark that the strategy proposed above is quite similar to the mod¬ 
ified equal time Wigner function formalisrrp^ used in [HI Sec. 5.2]. There, 
one subtracts from the two-point function not the Hadamard parametrix, 
but the vacuum two-point function, multiplied by a parallel transport. This 
may be appropriate in the 1-1-1 dimensional setting considered there, as the 
remainder is continuous (though not continuously differentiable) for equal 

^^Note, however, that the Hadamard coefficients Vk that appear on the r.h.s. of 0 
are recursively defined by an inhomogeneous ODE, cf. |20| for details. 

^"^Furthermore, in the context of the back-reaction problem of the semi-classical Einstein 
equation, a similar approach was taken in [33] for the evaluation of the stress-energy tensor. 
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times. However, it easily follows from the results in m. that in 3+1 dimen¬ 
sions the remainder would not be continuous. In particular, the logarithmic 
divergence that is responsible for the charge renormalization ambiguity is 
not canceled, so one expects the problems with greatly enhanced values at 
intermediate times to persist. 

5 Conclusion 

To summarize, we have introduced a method to compute the expectation 
value of the current density 0 in QED in external time-dependent poten¬ 
tials at intermediate times, based on a point-splitting renormalization. We 
showed that it is equivalent to the expression for the current density derived 
in the QKE framework [18]. We also corrected an expression due to Berber 
for the linearization of the vacuum current density in the external potential. 
Our non-perturbative results for the Sauter pulse Q are compatible both 
with the corrected version Q of Berber’s general expression in the case of a 
sub-critical peak field strength and with the asymptotic results of Narozhnyi 
and Nikishov [M] . 

We proposed a natural explanation of the enhanced quasi-particle densi¬ 
ties at intermediate times for slowly varying sub-critical homogeneous elec¬ 
tric fields [5] . The argument used the fact that the current is completely fixed 
up to the charge renormalization ambiguity. We also sketched a framework 
in which no difference of scales of the dynamical quantity at intermediate 
and asymptotic times is to be expected, and which is in principle applicable 
to generic external potentials. The development of this scheme is a topic for 
future work. 

Acknowledgements 

I would like to thank Bergey Gavrilov, Holger Gies, Jan Bchlemmer and 
Half Bchiitzhold for helpful discussions or hints to the literature. Parts of 
this work were done at the Faculty for Physics of the University of Vienna, 
supported by the Austrian Bcience Fund (EWE) under the contract P24713. 

A Comparison of point-splitting and QKE 

To fix conventions, we choose signature (+, —, —, —) and the standard Dirac 
matrices. The covariant derivative is defined as 

— dfj, + 

We construct a Klein-Gordon type operator 
P = — m) = 


13 


As only one independent component of F^i, is non-vanishing, namely Fq^, 
we choose two linearly independent eigenvectors 


= 


/o\ 
1 
0 

v-iy 


= 


0 

1 

Vo/ 


of 7 ^ 7 ^ with eigenvalue 1. On these, the operator P acts as the scalar 
operator 

P = + ieFos + 

We assume = h 3 ^A(x'^) and look for solutions gp of the form gp{x) = 
fp{x^)e^^F. Hence, fp fulfills 

fp++ ieE) fp = 0, (12) 

with TT as in Q. Assume that A{t) is switched on at some finite time and 
let fp{x^) = in the past. Properly normalized spinorial solutions to 

the Dirac equation are then given by 

Tp,S 




{-if - m)5fpr'* 


Using that 

one easily obtains 




j;r*f* = 7°-73, 


(^i - ^Dl/pl^ - \ fp\ - 


{2Trfcjp{u!p -pa) 


(13) 


On the other hand, for the vacuum modes evaluated at the shifted momen¬ 
tum, we have 


(14) 


Furthermore, it follows from (12) that 

dt\f/ = 2?fi{fpfp), 

dtifpfp) = \fpf - {^^l + ieE)\fpf, 
dt\f/ = -2ujliR{fpJp) - 2eEQ{fpJp). 
Comparison of ([^ with ( [T^ , (14) suggest that 

vrs {^l-'^i)\fp\'^-\fpf-‘^'^s’^ifpfp) 


(15) 


'^3 ■ 

—n H—=n =- 

lOtt eh 


2ojp{ujp - pfj 


(16) 
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Indeed, using (15), one can show that with 


n = 1 + 


TTslw^-l/pl + ^\fp\ ) + 2uj^Q{fpfp) 


2uip{uip P 3 ) 


( 17 ) 


one fulfills not only (16) but also the system of equations 0 , with the 
correct initial conditions. Furthermore, it is straightforward to verify that 


(27r)3 y 127r2 


+ 0 ( 2 ). (18) 


This proves the equivalence of ([^ and ([^. 


B The correction of Serber’s result 


Let us recall the arguments of Berber [2^ in his derivation of the vacuum 
current at linear order in the external field. He finds that, in units where 
m = 1, 




a 


r^/2 


327r5 


cos^ / A{x — x',ilj)ODJ^[x')d‘^x', 


wher 


/ 


A{x,ijj) = / e log(l — cos'!/;)A: ^d^fc. 


(19) 


cf. equations {13} and {14} (equation numbers in curly brackets refer to 
[26]). It is then shown that A vanishes outside of the light cone and is given 
by 

27r3 cos V’ 


A(a;, 'll)) = f j^(2A:-\/x 2/cos V’)A: ^dfe 

V x2 Jl 


( 20 ) 


inside, cf. equation {17}. To determine the value on the boundary of the 
forward light cone, assumed to be of the form /(r, ip)[6{t—r)-\-6{t+r)], Berber 
proceeds as follows: He first argues that the integral of (poj) over x^, for any 


fixed spatial x*, vanishes. He then computes the same integral for (19), finds 


a non-vanishing result, and concludes that it must be due to a non-vanishing 
/, i.e., a contribution localized on the boundary of the light cone. However, 
both these integrals are not computed correctly. Regarding the vanishing of 
the time integral over (20), Berber first introduces the variable s = 


and correctly notes that |35l Beet. 13.6] 

poo 

/ Ji{2ks/ cos'i/’)/'\/ + r^ds = Ii{k/ cos'4>)Ki(k/ cosi/j), 

Jo 2 2 

^®To obtain the retarded response, one has to replace in the log and the denominator 
by k^+ieko. Berber does not do this at this stage and obtains the mean of the retarded and 
advanced response, only selecting the retarded response at the very end (and multiplying 
it by 2). This explains the missing factor of 2 in (201 with respect to ([^. 
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but then argues that this vanishes, due to Ki (z) = 0. However, it is well- 

_2 

known [36l Sect. 10.2] that Ki{z) = y^vr/ {2z)e~^. Using computer algebra, 
the time integral of ( 20 ) can be computed explicitly, yielding 




=- jcos^?/)-!-(grcosV'- iCos^V')e ^UcosV' 


-I- Ei(—2r/ cos V’) 


where Ei denotes the exponential integral. 
Serber claims that 


For the time integral of (19), 




87r2 

.3 . 


/ sin(fcr) log(l-|- 7 /c^ cos^'!/^)/A:^dA: 

r Jo 


2tt'^ cos^ ijj 


i: 


^—2krl cos 


^/k^dk, 


the step from the first to the second line being the result of a deforma¬ 
tion of the integration path to the imaginary axis. However, this is not 
straightforward, as the closing of the path at infinity diverges, due to the 
sin. Numerical integration shows that integrals in the first and the second 
line do not coincide, so the path deformation seems to have been performed 
incorrectly. Performing the integral in the first line by computer algebra, one 
finds that = /([ 2 o|i(u V')- Hence, no supplementary term localized 

on the boun^ry of the li^t cone is present. 
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